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We discuss the dynamics of ultracold atoms in an optical 
potential accelerated by gravity. The positions and widths 
of the Wannier-Stark ladder of resonances are obtained as 
metastable states. The metastable Wannier-Bloch states os- 
cillate in a single band with the Bloch period. The width of 
the resonance gives the rate transition to the continuum. 



I. INTRODUCTION 

Optical lattices are a cornerstone of much of atom op- 
tics. They have been successfully used to form matter- 
wave diffraction gratings, were at the heart of the discov- 
ery of subrecoil cooling, and were manipulated to study 
atomic Bloch oscillations ||l|, and dynamic instabilities 
and chaos. ||^ They are also used as an outcoupling mech- 
anism in a mode-locked atom laser S , and are central to 
several quantum computing schemes using neutral atoms. 
Q As atom optics rapidly moves toward applications, fre- 
quently listed possibilities include gravimeters and grav- 
ity gradiometers. Gravity is also an important and nat- 
ural way to accelerate atoms, so as to coherently reduce 
their de Broglie wavelength for nanofabrication applica- 
tions. Hence, there is considerable interest in studying 
the combined effects of an optical lattice and gravitation 
on the dynamics of ultracold atoms. 

The study of particles in periodic potentials has of 
course been a major tenet of condensed matter physics, 
and it will therefore come as no great surprise that our 
analysis relies heavily on the use of its tools. In partic- 
ular, delocalized Bloch states and band structure anal- 
yses are known to be ideal for the study of electrons in 
crystals. For deep potentials, though, the use of localized 
Wannier states is often preferable and indeed, has proven 
useful e.g. in the study of the dipole-dipole interaction 
in optical lattices. |g] When a linear potential is superim- 
posed to the underlying periodic structure, the situation 
is less clear: the localized Wannier states become delocal- 
ized, and there is some debate as to which representation 
is then most useful. One possibility is to use a so-called 
Wannier-Bloch states representation, with an emphasis 
on the delocalized character of the atomic wave func- 
tions. Alternatively, the Wannier-Stark approach em- 
phasizes the localized aspect of the wave functions. A 
third alternative consists in treating the periodic poten- 
tial as a perturbation on top of the gravitational field. 



Clearly, each approach has its pros and cons, depending 
upon the depth of the lattice compared to the gravita- 
tional potential. For deep lattices, though, the Wannier- 
Stark analysis seems preferable. In this approach, the 
states of the atoms form a sequence of resonances sepa- 
rated by equal energy intervals referred to as Wannier- 
Stark ladders. |^-^ Such ladders were recently observed 
in accelerated optical lattices. PJlO) 

In this paper we calculate the Wannier-Stark ladders 
induced by gravity on atoms in optical lattices, using a 
theoretical approach recently developed in Ref. ||ll[ and 
based on a Floquet representation of the Wannier states. 
This technique, which leads to a well-defined eigenvalue 
problem for the determination of the Wannier-Stark lad- 
der, was developed for the more general situation of time- 
periodic potentials, but is also perfectly adapted to time- 
independent potentials. 

Section II introduces our physical model, and briefiy 
reviews the definitions of the Bloch and Wannier states, 
which are the eigenstates of the system in the absence of 
a gravitational field. Section III outlines the Floquet ap- 
proach to the determination of the Wannier-Stark ladder, 
discussing the onset of a ladder of complex "resonances" 
associated with metastable localized states of the atoms. 
Selected numerical results are presented in section IV. 
Finally section V is a summary and outlook. 



II. PHYSICAL MODEL 

We consider the situation of ultracold atoms of mass 
M trapped in the optical potential resulting from two far- 
off resonant and counterpropagating optical fields of fre- 
quency ojl and wave number k^ — 2-k /\l. The atoms are 
assumed to be two-level atoms, and we neglect all compli- 
cations resulting from the possible existence of multiple 
ground states and the associated optical pumping. After 
adiabatic elimination of their upper electronic states, the 
atoms are therefore described as scalar particles subject 
to the periodic potential 
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where Vi is the Rabi frequency of the atomic transition 
of frequency loa under consideration. In the presence of 



gravity, the atomic Hamiltonian of the system is given 

by 



H ^ Ho + Mgz, 
where Mg is the force of gravity and 

^o^l^ + m 



(3) 



(4) 



We consider a one-dimensional situation where the lasers 
generating the optical lattice are along the z-axis, so that 
[z,p] =ih. 

The band structure associated with the Hamiltonian 
Hq is of course well-known. It is displayed in Fig. 1 for 
reference. It shows the energy-wave number dispersion 
relation £'„(fc) in the first Brillouin zone [|2|, n being 
the band index. The associated energy eigenstates of Hq 
are the familiar Bloch states 

Ho4)nM^) = En{k)(l)n,k{z). (5) 

They are conveniently expressed as 

(j)n,k{z)=eyi'p{ikz)xn,k{z), (6) 

where the functions Xn.k{z) are periodic in z. They sat- 
isfy the eigenvalue equation 



Ho''\„,k{z) = £'„(fc)x„,fe(z), 
the Hamiltonian TJq being given by 



H, 



(.) ^ {p + hkr 
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From Eq. (M), the evolution of the Bloch wave func- 
tions Xn,fc(2) is simply 



Xn,k{z,t) = exp[-i£;„(fc)i/?i]xn,fc(2,0). 



(9) 



It is furthermore possible to take advantage of the pe- 
riodicity of the optical potential to expand Xn,k{z) on a 
lattice basis {|£)}, € integer, as the Fourier series 



„,fe(z) = \/ — V 4"^'^ exp(2i£fciz), (10) 



where we have used the coordinate representation of 



{z\(.) 



exp(2z£fcLz), 



(11) 



and we recognize the primitive reciprocal vector h = 2kL 
in the exponent. 

In prevision of our subsequent discussion, we introduce 
at this point the characteristic time 



Tb = 



2hkT 



Mg ' 



(12) 



which will later on be identified as the "Bloch period," 
as well as the associated "Bloch phase" 



e„(fc) - E,,{k)TB/n. 



(13) 



Since the energies En{k) always appear in the form of a 
phase factor in the evolution of the system, it is in many 
cases convenient to plot them in a reduced phase inter- 
val. This representation will be useful in identifying the 
resonances appearing in the presence of a gravitational 
potential. 

We conclude this section by recalling that the Bloch 
states are completely delocalized, with the atomic wave 
functions extending over the whole lattice. For deep lat- 
tices, it is oftentimes more convenient to use instead lo- 
calized Wannier states. They are constructed as linear 
superpositions of the Bloch states 4>n,k{z) as 



y:>n,m{z) = / exp{-iRmk)(j)n,k{z)dk, (14) 

where Rm is the position of the m-th lattice site. 



R„ 



m-K/kL, 



(15) 



m is an integer, and we recognize that Tr/fc^ is just the 
primitive lattice period a of the optical potential. As 
is well known, the Wannier function ipn,m{z) describes 
an atomic wave packet localized at the m-th site of the 
optical lattice. 



III. WANNIER-BLOCH STATES AND 
WANNIER-STARK STATES 

We now turn to a discussion of the impact of grav- 
ity on the atomic dynamics in the lattice. We proceed 
by eliminating the gravitational contribution from the 
Hamiltonian (|3|) via the unitary transformation 



il){z,t) = eyi\i[—iMgzt/h]tl){z^t). 



(16) 



The evolution of ^{z, t) is governed by the time- 
dependent transformed Hamiltonian 



nt) 



{P - Mgtf 
2M 



V{z), 



(17) 



which is nothing but the Hamiltonian Hq with momenta 
shifted by Mgt. From this observation, it follows that its 
time evolution is governed by the equation 

n{t + T) = exp{iMgTz/n)H{t) exp{-iMgTz/n). (18) 

Eq. (^^ provides an immediate understanding of Bloch 
oscillations: After a time Tb — 2hkL/Mg, the momen- 
tum shift appearing in the Hamiltonian dl8) is precisely 
equal to 2kL, the primitive reciprocal lattice vector. An 
atom starting at the edge of the first Brillouin zone is 



therefore uniformly accelerated into the next Brillouin 
zone, or equivalently undergoes an umklapp process in 
the first zone. In optics terms, this means that the Bragg 
scattering condition is fulfilled, with scattering of the 
atom from the state of momentum k^ to —k^. This jus- 
tifies identifying Tb as the Bloch period as we have done 
in section II. 

With this intuition at hand, we introduce the Wannier- 
Bloch states ipn,k{z) as those states which are strobo- 
scopic eigenstates of the time-dependent Hamiltonian 
( p^ ) at successive intervals of Tb, that is, that satisfy 
the eigenvalue equation 



U(TB)'4'n,kiz) = exp 



-^En{k)TB 



with 



U{Tb) =cxp -- / Hdt 



V'«,fe(2). (19) 



(20) 



Despite the fact that we use the same notation En{k) 
for simplicity, it is important to realize that the quantities 
appearing in Eq. (n9h are not the same as the band ener- 
gies of section II. Indeed, due to the fact that the Hamil- 
tonian in the presence of gravity is time-dependent, the 
En{k) are only quasi-energies, or "resonances." We will 
see shortly that they are complex, their imaginary part 
leading to a finite lifetime, and form a series of ladders 
of equidistant levels. 

The eigenstates ipn,k{z) present the considerable ad- 
vantage that the only momentum involved at the end 
of their stroboscopic evolution is the primitive reciprocal 
lattice vector 2kL- That is, U{Tb) commutes with the 
momentum shift operator exp[— iTrp/Tifci], and hence it 
is still possible to use a Bloch-like quasi-momentum rep- 
resentation to stroboscopically diagonalize the problem. 
Rewriting then, in analogy with Eq. M, the Wannier- 
Bloch functions as 



tpn,k{z) = eyip{ikz)x„Mz), 



(21) 



where Xn,k{z) = Xn.kjz + 7r/fci), leads to the Bloch-like 
eigenvalue equation pT| 



i^^'hnA^) 



exp 



-E^{k)TB 



Xu,k{z), (22) 



where the operator U^ ' is given by 

and 

U^'^HTb)^ 



(23) 



Texp 



Tb 



{p + hk- Mgtf 



hJo\ 2M 

where T stands for time-ordered. 



+ V{z) dt 



(24) 



From the analysis so far, it woul d ap pear that we need 
to solve the eigenvalue problem (22) for a large num- 
ber of initial times between and Tb to obtain the full 
system evolution. This, however, is not the case, as 
the quasi-energy spectrum En{k) is completely degen- 
erate for the time-independent potential considered here, 
En{k) — En- This can be seen by considering the evolu- 
tion of a Wannier-Bloch state for a fraction T = Tb/s of 
the Bloch period, where s is an integer. In that interval 
the operator 



Ut 



-liki^z I S 



U{T), 



with 



Ut=U, 



(25) 



(26) 



transforms the Wannier-Bloch state ipn.k{z) into another 
Wannicr-Bloch state of quasi-momentum k' = k — Ikj^j s, 

UTlpn,k{z)^^n,k'{z). (27) 

Combining Eqs. (^ and ( p7| ) yields readily 

Ulpn,k'{z) = iiiiTlpn,k{z) = Urii'ipnMz) 

= Ut exp 



'TEn{k)TB 

n 



exp 



-J.En{k)TB 



i^n.kiz) 
1pn,k'iz) 



(28) 



where we have used Eq. ( |26| ) to commute U and Ut- 
This completes the proof that for the time-independent 
potential V{z) at hand, the quasi-energy spectrum is de- 
generate. 

As a result, the continuous time evolution of the 
Wannicr-Bloch states is simply given by 

^n,k{z,t) ^ exp{-iEnt/Ti)ipn^k-Mgt/h{z,0). (29) 

The calculation of the quasi-energies £'„ follows the 
procedure discussed in Rcf. ||ll| . It consists in finding the 
metastable Wannier-Bloch states by solving the eigen- 
value problem d23) assuming a truncation of the momen- 
tum space. The representation of the operator U^'^' in 
this truncated basis is nonunitary, and for this reason 
the "energies" En are complex. 



En -^ En— i^nl'2- 



(30) 



Their real part give the positions of the "resonances" , 
while their imaginary part give their widths, or lifetimes. 
This is to be contrasted to the gravitation-free situation, 
where the eigenvalue problem is tridiagonal in momen- 
tum space, and the eigenvalues are then real. Eq. ( [2^ ) 
shows clearly how the acceleration of gravity removes the 
discrete, tridiagonal nature of the problem and makes it 
continuous in momentum space. 

We have previously discussed how the Bloch func- 
tions approach puts the emphasis on dclocalized states. 



Clearly, this same characteristic also holds for the 
Wannier-Bloch states discussed so far. As in the 
gravitation-free situation, it is sometimes useful to con- 
sider states of a localized nature instead. In complete 
analogy with the Wannier functions, we therefore intro- 
duce the Wannier-Stark states as 



^n^z) 



e-Kp{-iRmk)ip„^k{z)dk, (31) 



-fct 



where as we recall Rm = rmr/k]^ is the position of the 
77i-th lattice site. Taking into account the temporal evo- 
lution (E9h of the Wannier-Bloch states gives the contin- 
uous time evolution of the Wannier-Stark wave functions 



*n,m(2, t) = exp[--{E„ + mT:Mg/kL)t]^n,m{z)- (32) 

h 

This shows that the Wannier-Stark states constitute a 
ladder of quasi-stationary functions characterized by sets 
of discrete levels of energies 



En,7n ^ En+ imrMg/kL, 



(33) 



of equal spacings irMg/k]^. These energy levels are illus- 
trated in Fig. 2. 



IV. RESONANCES 

In this section, we present selected results from the nu- 
merical study of gravitationally-induced Wannier-Stark 
resonances for the case of shallow optical lattices. The 
calculations proceed by solving the eigenvalue equation 



r„(fc) 



exp 



2^ rs)exp(--^„(fc)rB)x«,fc(^), (34) 



where W'*-') is the operator (p3), truncated in momentum 
space as previously discussed. 

To diagonalize the matrices VV^''^ (£',£), we first evalu- 
ate the operator U'^^^Tb) using the formula 



N-l 

U^^){TB) = r\{ exp 



'*' ■H^''\t)dt 



where 



n^'^Ht) 



{p + hk- Mgtf 
2M 



V{z), 



(35) 



(36) 



and t.j+i - tj = At = Tb/N {N » 1). The exponen- 
tial factors in ( p5| ) are approximated by a Dyson series 
truncated at second order. 



exp 



n'^^\t)dt 



! + (--) / n^^\t)dt 



+ (--)' / H^''Ht)dt / 'H^^\t')dt'. (37) 

The analytic evaluation of the matrix elements of the 
Dyson series (^) is a straightforward procedure. We 
just give the result of the second term: 



(£'| f '^' H^''Ht)dt\£) ^ 
Jtj 

E„T, 



EbTb 

N 



21- 



2j + l 

N 



3iV3 



6i> 



VqTb 

AN 



{Se\e+i + Si'^e-i) , 



(38) 



where Er — TT?k\/2M is the recoil energy, k — k/k^ and 



we have used the definition (h2|) of the Bloch period. For 
fixed K, Eq. ( |37| ) gives the elements of a symmetric ma- 
trix whose truncation causes no problem. We proceed by 
first truncating it, and by an iterative process of succes- 
sive time-ordered intervals we then evaluate the evolution 
operator W'^^Tb) using Eq. (pSf). The matrix elements 
of the operator e^'^^^^^ are given by 



{e\e 



-2iki^z 



1^) = -5^/- 



(39) 



whose truncation generates a nonunitary operator. Tak- 
ing the dimension of the matrices to be 2R -I- 1, we have 
found that for sufficiently large i?, typically, 20 to 40, 
the position of the resonances becomes independent of 
the dimension of the truncated momentum space. 

The simulations were carried out for the case of ^^ Rh 
atoms |3J trapped in an optical lattice generated by a far- 
off resonance optical standing wave of wavelength A = 
850nr7T,. The depth of the lattice was taken to be Vb = 
—2.1 En, where En « 3.177 kHz is the recoil energy. This 
choice of parameters corresponds to a level separation of 
AE^:s = TrMg/kL = 0.2854485£;j^ in the Wannier-Stark 
ladder, see Eq. ([3^). 

Fig. 3 shows two resonances, folded in the interval 
— TT < 6,1 (/c) < n where Qn{k) = En{k)TB/h (modulo 
27r), see Eq. ([l3|). They are degenerate in k, as dis- 
cussed in section III. (Recall that only degenerate solu- 
tions are associated with resonances.) The real part of 
these quasienergies are Eq = 0.10771571£^/{ and Ei — 
1.024990i?ii;, respectively, corresponding to the Bloch 
phases Qq = 2.37112 (solid line) and Oi = -2.57103 
(dotted line). The widths of the resonances, given by the 
imaginary part of the degenerate eigenvalues of Eq. (|34|), 
are shown in Fig. 4. For the relatively shallow lattice 
considered here, the highest resonance is probably too 
wide, with ri/2 = 7.182 x lO^^^;^^ (dotted line), to be of 
practical use. In contrast, the width of the most stable 
resonance, ro/2 = 2.17378 x IO'^^Er (sohd line), is three 
orders of magnitude narrower than the Wannier-Stark 
ladder separation AE^^s- As further discussed in the fol- 
lowing section, such narrow lines might be of consider- 
able interest in the development of atom optical gravity 
gradiometers. 



Despite the long lifetime of the metastable Wannier- 
Stark states, the fact that they are finite indicates that 
atoms initially localized in a lattice site eventually tun- 
nel to neighboring wells. This is shown in Fig. 5, which 
plots on a logarithmic scale the absolute squared expan- 
sion coefficients |c^"' ~ p of the Bloch states, see Eq. 
(O), for the lowest band n = 0. The solid line corre- 
spond to the gravity-free situation, while the dotted line 
correspond to gravity. When the gravity force is applied 
the metastable states tunnel through the optical poten- 
tial hills in the downward direction (negative z-axis), so 
the contribution of negative 2£kL vectors increases. In- 
deed, there is a general shift of all the coefficients c^"' 
due to the action of the gravity field. 

We have also evaluated the position and width of the 
resonances as a function of the well depth Vb of the opti- 
cal lattice. Figs. 6 and 7 show the position and width of 
the lowest resonance, Eo and ro/2 respectively, in units 
of Eh. As we can see in these plots, the position Eq goes 
up as the well depth decreases, while, the width ro/2 
decreases very rapidly for highest potential depths. Fi- 
nally, Figure 8 is a plot of the energy gap Egap between 
the lowest resonance and continuum states. 



V. SUMMARY AND CONCLUSIONS 

For a number of potential applications of atom optics, 
such as e.g. nanolithography and atom holography, it is 
desirable to accelerate the atoms so as to produce coher- 
ent beams of shorter de Broglie wavelength. This can 
be achieved either via gravity, or in lattices with ramped 
detuning between the two counterpropagating fields, as 
already mentioned. Since eventually, the atoms will then 
escape the trap, matter- wave acceleration is also of much 
relevance in the study of output couplers for atom lasers, 
as was realized for instance the mode-locked laser of Ref. 
0. In addition, atom optics with accelerated beams of- 
fers novel approaches to interferometry. For instance, us- 
ing optical fields to hold atoms against gravity may lead 
to the realization of accelerometers of unprecedented sen- 
sitivity. 

These new developments require a detailed under- 
standing of the combined effects of light and gravity on 
the motion of ultracold atoms. Using a recently devel- 
oped Floquet-type analysis |I^ that removes many of 
the question marks previously attached to the computa- 
tion of the Wannier-Stark ladder, we have shown that 
a detailed quantitative calculation of the position and 
width of these levels can now easily be obtained numer- 
ically. We illustrated how deep optical lattices lead to 
the existence of families of long-lived, equidistant quasi- 
energy levels that should be of much interest in grav- 
itation gradiometers and accelerometers. As the wells 
become shallower, however, the number of quasi-bound 
levels decreases, and their lifetime shortens. 

We note that the applications we have in mind are 



likely to eventually require the use of atomic samples, 
most likely Bose-Einstein condensates, for which atom- 
atom interactions and the associated nonlinear phase 
shifts can no longer be ignored. This is likely to influence 
the Wannier-Stark spectra in interesting ways, that we 
plan on studying in the future. 
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FIG. 1. Energy bands En(k) (in units of Er) of an opti- 
cal lattice as function of the quasimomentum k in the first 
Brilouin zone [— fci,,fci,l for V{) — —2.1Er. 
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FIG. 2. Tilted periodic potential and Wannier-Stark lad- 
ders, in units of the recoil energy Er, as functions of z for 
Vo = -2.1Er. 
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FIG. 3. Numerical results for resonances folded in the in- 
terval -TT < e(fc) < TT with e(fc) = E„{k)TB/h (mod- 
ulo 2n) for Vo = -2.1Er. Go = 2.37112 (solid line) and 
01 = —2.57103 (dotted line). They correspond to energies 
Eo = 0.10771571£;_R and Ei = 1.024990^^, respectively. 



FIG. 4. Numerical results for the width of the resonances 
r„(fe)/2 (in units of Er). ro/2 = 2.17378 x 10~*Er (solid 
line) and ri/2 = 7.182 x W'^Er (dotted line). 
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FIG. 5. Absolute squared coefficients |c^"' p on a loga- 
rithmic scale for band n = and quasimomentum A; = in 
the field-free case (solid line) and gravity case(dotted line). 
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FIG. 6. Lowest resonance _Eo (in units of Efi) as function 
of the potential depth Vo . 
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FIG. 7. Width ro/2 (in units of Er) as function of the 
potential depth Vo- 
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FIG. 8. Energy gap Egap (in units of Eb) between the low- 
est resonance and continuum states as function of the poten- 
tial depth Vq. 



